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ABSTRACT 


The  present  paper  continues  the  study  of  balanced  treatment  incomplete 
block  (3TIB)  designs  initiated  in  [l]-[4].  This  class  of  designs  was 
proposed  for  the  problem  of  comparing  simultaneously  p  >  2  test  treat¬ 
ments  with  a  control  treatment  when  the  observations  are  taken  in  blocks 
of  common  size  k  <  p+1.  The  conjectured  minimal  complete  class  of 
generator  designs,  a  catalog  of  admissible  designs,  and  tables  of  optimal 
designs  are  given  for  p  =  6,  k  =  3.  The  efficiency  of  the  optimal  BTIB 
design  relative  to  that  of  replications  of  a  BIB  design  is  computed  for 
situations  in  which  both  provide  the  same  probability  guarantee  for  the 
multiple  comparisons  with  a  control  problem. 

Key  words  and  phrases:  Multiple  comparisons  with  a  control,  balanced 
treatment  incomplete  block  (BTIB)  designs,  BIB  designs,  admissible  designs, 
S  -inadmissible  designs,  C -inadmissible  designs,  minimal  complete  class 
of  generator  designs,  optimal  designs. 


1. _ INTRODUCTION 


The  present  paper  continues  the  study  of  balanced  treatment  incomplete 
block  (BTIB)  designs  discussed  in  tl]-[4].  This  class  of  designs  was  pro¬ 
posed  for  the  problem  of  comparing  simultaneously  p  >  2  test  treatments 
with  a  control  treatment  when  the  observations  are  taken  in  blocks  of 
common  size  k  <  p+1. 

In  [1]  a  general  theory  of  BTIB  designs  was  developed;  in  [2]  optimal 
designs  were  given  for  the  cases  p  =  2,  k  =  2(1)6  and  p  =  3,  k  =  3. 

In  [3]  optimal  designs  were  provided  for  the  cases  p  =  4,  k  =  3  and 
p  =  5,  k  =  3  while  in  [4]  optimal  designs  were  given  for  the  case  p  =  4, 
k  =  4.  In  the  present  paper  we  give  optimal  designs  for  the  case  p  =  6, 

k  =  3;  these  optimal  designs  are  subject  to  the  same  qualification  as  those 

given  in  [3]  and  [4] — namely  that  they  are  optimal  relative  to  the  generator 
designs  known  to  us. 

The  reader  is  referred  to  [3]  for  the  definitions  of  inadmissibility, 

S-inadmissibility  and  C- inadmissibility  used  in  this  paper.  The  reader 

is  also  referred  specifically  to  Sections  2  and  3  of  [2]  and  Sections  1  and 

2  of  [3]  for  an  exact  statement  of  the  multiple  comparison  problem  under 

consideration,  expressions  for  the  BLUE's  of  the  treatment  effect  differences 

2  2 

aQ  -  cu  (1  <  i  <  p),  their  variances  (t  c  )  and  correlations  (p),  and 

an  expression  for  the  confidence  coefficient  (P)  associated  with  joint  one¬ 

sided  confidence  interval  estimates  of  the  aQ-ou  (1  <  i  <  p), 

2.  RESULTS  FOR  p  =  6,  k  =  3 

2.1  Conjectured  minimal  complete  class  of  generator  designs 

In  Table  2.1  we  list  for  p  =  6,  k  =  3  the  generator  designs  in  our 
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Table  2.1 

Conjectured  Minimal  Complete  Class  of  Generator  Designs  for  p  =  6,  k  =  3 
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conjectured  minimal  complete  class.  We  prove  in  the  Appendix  that  in 

that  class  is  C-inadmissible  for  b  >  22. 

2.2  Catalog  of  admissible  designs 

A  catalog  of  admissible  designs  has  been  prepared  based  on  the  conjectured 
minimal  complete  class  of  generator  designs  given  in  Table  2.1.  This  catalog 
is  given  in  Table  2.2  for  b  =  6,7,11,14,15,17,18,21,22,24(1)63. 

2.3  Tables  of  optimal  designs 

Optimal  designs  for  p  =  6,  k  =  3  are  given  in  Table  2.3  for  d/a  =  0.1(0. 1)1.0 
for  b  =  6-62.  Optimal  designs  that  achieve  a  specified  confidence  coefficient 
1-a  are  given  as  a  function  of  d/a  =  0.2(0. 2) 2.0  for  1-a  =  0.75,  0.80,  0.85, 

0.90,  0.95  and  0.99  in  Table  2.4.  These  designs  were  found  by  a  complete  computer 
search  among  all  admissible  designs. 

3.  EFFICIENCY  OF  THE  OPTIMAL  BTIB  DESIGN  RELATIVE  TO  THAT  OF 
REPLICATIONS  OF  THE  BIB  DESIGN  D2  FOR  THE  MCC  PROBLEM 

It  is  of  some  interest  to  compare  the  efficiency  of  the  optimal  BTIB  design 
relative  to  that  of  replications  of  the  BIB  design  D2  for  the  MCC  (multiple 
comparisons  with  a  control)  problem  when  each  is  used  to  achieve  the  same  con¬ 
fidence  coefficient  1-a  for  a  given  d/a.  We  have  computed  the  required  number 
of  blocks  for  the  BIB  design  and  given  the  results  in  Table  3.1  along  with  the 
corresponding  results  for  the  optimal  BTIB  design  (the  latter  being  abstracted 
from  Table  2.4);  also  given  in  the  table  is  the  ratio  of  the  number  of  blocks 
required  by  the  BTIB  design  to  that  of  the  number  of  blocks  required  by  the 
BIB  design.  This  ratio,  which  we  term  the  relative  efficiency,  must  be  less 
than  or  equal  to  unity  since  a  BIB  design  is  also  a  BTIB  design.  Thus  the 
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Table  2.2 

Catalog  of  admissible  designs— ^  for  p  =  6,  k  =  3 


No. 

of 

blocks 

(b) 

D1 

bl  =  6 

A(1)  =  2 
u  2 

=  0 

°2 

b2  =  7 

a(2)  =  1 
Ao  i 

a(2)  =  1 
A1  A 

D3 

b3  =  11 

A(3)  =  3 

Ao  i 

a<3)  =  1 

D4 

b4  =  15 

x(4)  =  5 

Ao  5 

44)  = 1 

°5 

b  =  10 

4s  1  ■  o 

45)  =  2 

xo 

i 

2 

T 

P 

6 

1 

0 

0 

0 

0 

E 

0 

1.5000 

0.000 

7 

0 

1 

0 

0 

0 

B 

0.8571 

0.500 

11 

0 

0 

1 

0 

0 

3 

0.4444 

0.250 

14 

0 

2 

0 

0 

0 

2 

0.4286 

0.500 

15 

0 

0 

0 

1 

0 

5 

E 

0.3273 

0.167 

17 

0 

1 

0 

0 

1 

B 

3 

0.6316 

0.750 

18 

0 

1 

1 

0 

0 

1 

2 

0.2813 

0.333 

21 

0 

0 

1 

0 

1 

3 

3 

0.2857 

0.500 

1 

0 

0 

1 

0 

1 

0.2637 

0.125 

22 

0 

1 

0 

1 

0 

2 

0.2222 

0.250 

24 

0 

2 

0 

0 

1 

E 

4 

0.3462 

0.667 

25 

0 

2 

1 

0 

0 

3 

0.2087 

0.375 

26 

0 

0 

1 

1 

0 

8 

2 

0.1875 

0.200 

27 

0 

1 

0 

0 

2 

0.5806 

0.833 

28 

0 

1 

1 

0 

1 

E 

0.2143 

0.500 

29 

0 

1 

2 

0 

0 

0.1714 

0.300 

30 

0 

0 

0 

2 

0 

E 

0.1636 

0.167 

— ■  For  each  number  of  blocks,  the  number  under  Cl  <  i  ^  5)  in  the 

body  of  the  table  is  the  frequency  f.  with  which  D.  appears  in  the 
5  1  1 

design  D  =  u  f .D. . 

i=l  1  1 
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Table  2.2  (continued) 


Mo. 

of 

blocks 

(b) 

°1 

b1  =  6 

A(1)  =  2 
u  1 

a[1)  =  0 

D2 

b2  =  7 

x(2)  =  1 
0  1 

A^2)  =  1 

°3 

b3  =  n 

x(3)  =  3 

a<3)  =  1 

D 

4 

b4  =  15 

xo  5 

44)  =  1 

D5 

b,  =  10 

<»  ■  • 
A<5)  =  2 

xo 

2 

T 

P 

31 

0 

3 

0 

0 

1 

3 

5 

0.2424 

0.625 

32 

0 

1 

0 

1 

1 

6 

4 

0.1667 

0.400 

33 

0 

1 

1 

1 

0 

9 

3 

0 . 1481 

0.250 

34 

0 

2 

0 

0 

2 

2 

6 

0.3158 

0.750 

35 

0 

2 

1 

0 

1 

5 

5 

0.1714 

0.500 

36 

0 

0 

1 

1 

1 

8 

4 

0.1406 

0.333 

37 

0 

1 

0 

0 

3 

7 

0.5581 

0.875 

0 

1 

0 

2 

0 

[1] 

3 

0.1317 

0.214 

38 

0 

1 

1 

0 

6 

0.1875 

0.600 

39 

0 

1 

2 

0 

1 

5 

0.1390 

0.417 

40 

0 

2 

1 

1 

0 

!m 

4 

0.1235 

0.286 

41 

0 

3 

0 

0 

2 

3 

7 

0.2222 

0.700 

0 

0 

1 

2 

0 

13 

3 

0.1191 

0.188 

42 

0 

1 

0 

1 

2 

6 

6 

0.1429 

0.500 

43 

0 

1 

1 

1 

1 

9 

5 

0.1197 

0.357 

44 

0 

2 

0 

0 

3 

2 

8 

0.3000 

0.800 

0 

1 

2 

1 

0 

12 

4 

0.1111 

0.250 

45 

0 

2 

1 

0 

2 

5 

7 

0.1532 

0.583 

0 

0 

0 

3 

0 

15 

3 

0.1091 

0.167 

46 

0 

2 

2 

0 

1 

8 

6 

0.1193 

0.429 

47 

0 

1 

0 

o 

4 

f 

9 

0.5455 

0.900 

0 

1 

0 

2 

1 

1 

5 

0 . 1064 

0.313 

48 

0 

1 

1 

0 

3 

3 

0.1731 

0.667 

0 

1 

1 

2 

0 

B 

4 

0.1015 

0.222 
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Table  2.3  (continued) 
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Table  2.3  (continued) 
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Table  2.3  (continued) 
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5406  0.8220  0.9145  0.9644  0.9871  0.9959  0.9989 
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Table  2.4 

Optimal  Design—^  to  Achieve  a  Specified  Confidence  Coefficient 

as  a  Function  of  d/a 
for  p  =  6,  k  =  3 


Confidence 

Coefficient 

(1-a) 

0.2 

0.4 

0.6 

0.8 

m 

ra 

D 

D 

1.8 

2.0 

b=1044 

b=261 

b=117 

b=66 

l 

B|jl 

m 

m 

m 

b=ll 

0.99 

0,1,2 

59,13 

0,1,4 

12,3 

0,1,5 

3,1 

0,1,4 

1,0 

m 

m 

1 

H 

yQH 

0,0,1 

0,0 

b=680 

b=171 

b=77 

b=44 

b=29 

b=22 

b=15 

b=ll 

b=ll 

b=ll 

0.95 

0,1,3 

36,10 

0,1,4 

6,3 

0,1,0 

4,1 

0,1,2 

1,0 

0,1,2 

0,0 

0,1,0 

1,0 

0,0,0 

1,0 

0,0,1 

0,0 

0,0,1 

0,0 

0,0,1 

0,0 

b=522 

b=131 

b=59 

b=33 

b=22 

b=15 

b=ll 

aa 

b=7 

0.90 
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1,0 
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0,0 

mi 

SI 

HW 
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0,1,0 

0,0 

b=428 

b=108 

b=48 

b=29 

1 

mm 

m 

b=7 

b=7 

0.85 

0,1,1 

20,11 

0,1,1 

4,3 

0,1,1 

2,0 

0,1,2 

0,0 
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a 

BBS 
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yiH 

0,1,0 

0,0 

0,1,0 

0,0 

b=360 

b=90 

b=40 

b=25 

b=18 

b=ll 

b=ll 

b=7 

b=7 

0.80 

0,1,3 

16,8 

0,1,3 

2,2 

0,2,1 

1,0 

0,2,1 

0,0 

0,1,1 

0,0 

0,0,1 

0,0 

0,0,1 

0,0 

0,1,0 

0,0 

Wm 

Bl 

0,1,0 

0,0 
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b=78 

b=36 

b=21 

b=14 

b=ll 
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0,1,4 
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0,1,1 
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1,1 

m 

1 
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0,0 

o,i,o 

0,0 

0,1,0 

0,0 

0,1,0 

0,0 

— ^The  "matrix"  in  each  cell  is 
5 


f  1  ’ f  2  * f  3 ' 


where  D  =  u  f.D.  with 


i=l 


i  l 


b  =  £  f.b,  is  the  optimal  design  for  the  given  value  of  1-a  and  d/a. 

i=l  1  1 


Table  3.1 


Efficiency  of  a  BIB  Design  Obtained  by  Replicating  D^ 
Relative  to  the  Optimal  BTIB  Design 
as  a  Function  of  d/o 
for  p  =  6 ,  k  =  3 


Confidence 

Coefficient 

(1-a) 

d/a 

0.2 

■B 

0.6 

0.8 

1.0 

0.99 

b  =  1044 
b  =  1253 
0.8332 

b  =  261 
b  =  315 
0.8286 

b  =  117 
b  =  140 
0.8357 

b  =  66 
b  =  84 
0.7857 

b  =  44 
b  =  56 
0.7857 

0.95 

b  =  680 
b  =  791 
0.8597 

b  =  171 
b  =  203 
0.8424 

b  =  77 
b  =  91 
0.8462 

b  =  44 
b  =  56 
0.7857 

b  =  29 
b  =  35 
0.8286 

0.90 

b  =  522 
b  =  588 
0.8878 

b  =  131 
b  =  147 
0.8912 

b  =  59 
b  =  70 
0.8429 

b  =  22 
b  =  28 
0.7857 

0.85 

b  =  428 
b  =  469 
0.9126 

b  =  108 
b  =  119 
0.9076 

b  =  48 
b  =  56 
0.8571 

b  =  29 
b  =  35 
0.8286 

b  =  18 
b  =  21 
0.8571 

0.80 

b  =  360 
b  =  385 
0.9351 

b  =  90 
b  =  98 
0.9184 

b  =  40 
b  =  49 
0.8163 

b  =  25 
b  =  28 
0.8929 

b  =  18 
b  =  21 
0.8571 

—  The  top  number  in  each  cell  is  the  smallest  number  of  blocks 
required  by  the  optimal  BTIB  design  to  achieve  1-a;  this  number 
is  taken  from  Table  2.4. 

2/ 

—  The  middle  number  in  each  cell  is  the  smallest  number  of  blocks 
required  by  the  BIB  design  obtained  by  replicating  D»  to  achieve 
1-a. 


3/ 

—  The  bottom  number  in  each  cell  is  the  relative  efficiency,  i.e., 
the  ratio  of  the  top  number  in  each  cell  to  the  middle  number  in 
that  cell. 
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optimal  BTIB  design  never  requires  more  blocks  for  the  MCC  problem  than  does 
replications  of  the  BIB  design;  the  smaller  the  ratio,  the  more  efficient 
is  the  optimal  BTIB  design. 

We  note  from  the  table  that  as  1-a  approaches  unity  for  fixed  d/a 

the  relative  efficiencies  decrease  (i.e.,  the  BTIB  design  becomes  relatively 

more  efficient)  while  as  d/a  approaches  zero  for  fixed  1-a  the  relative 

efficiencies  increase  (i.e.,  the  BTIB  design  becomes  relatively  less  efficient) 

This  apparent  inconsistency  can  be  explained  as  follows:  For  1-a  approaching 

unity  when  d/a  (d/a  >  0)  is  fixed,  it  is  known  that  (in  the  limit)  the 

2 

optimal  design  is  the  one  that  minimizes  x  ;  for  d/a  approaching  zero  when 
1-a  (1-a  <  1)  is  fixed,  it  is  known  that  (in  the  limit)  the  optimal  design 

is  the  one  that  maximizes  p .  These  two  limiting  optimal  BTIB  designs  are 
different,  and  the  limiting  relative  efficiencies  are  bounded  away  from  zero 
and  one,  respectively. 
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APPENDIX 


Proof  of  C-inadmissibility  of  for  b  >_  22 

We  note  that  D.  is  admissible  for  b  =  6,  and  D,  u  D  is  admissible 
l  14 

for  b  =  21.  We  further  note  that  D  u  D2  (b  =  13,  X  s  3,  \  *  1)  is 
S-inadmissible  wrt  02(11,3,1);  u  (17,5,1)  is  S- inadmissible  wrt 

D^  (15,5,1)  and  u  (16,2,2)  is  S-inadmissible  wrt  2D2  (14,2,2). 

Therefore  it  only  remains  to  show  that  every  design  of  the  form  f^D^  u  ^4^4 
is  inadmissible  except  (f.  ,f2)  =  (1,0)  or  (1,1)  which  are  admissible. 

Case  1  (f^  =  2m  with  m  ^  1)  :  We  shall  show  that  D  =  2mD^  u  f4D4  is 
inadmissible  wrt  D'  =  mD3  u  f^  for  m  >,  1,  f^  £  0.  Let  (b  ,XQ ,  A  ,0  ,t2) 
((&’ jAq.X^p' ,t  ))  denote  the  parameters  associated  with  D  (D1).  Then 
we  have  b  =  12m  +  15f4>  \Q  =  4m  +  Sf^,  X  =  f4  and  b'  =  11m  +  15f  , 

XQ  =  3m  +  5f4,  X^  =  m  +  f  .  Thus  b  >  b'.  We  have  p  <  p1  iff 
f 4/ ( 4m  +  6f4)  i  (m  +  f|+)/(4m  +  6f4)  which  holds  with  strict  inequality  when 
m  >  1.  Finally,  t2  ^  t’2  iff 

4m  +  6f4  4m  +  6f4 

(4m  +  5f4)(4m  +  llf4)  -  (3m  +  5f  )(9m  +  Ilf  ) 


which  holds  with  strict  inequality.  This  completes  the  proof  for  Case  1. 

Case  2  (f^  =  2m+l  with  n  >  1):  We  shall  show  that  D  =  (2m+l)D^  u  f4D4 
is  inadmissible  wrt  D'  =  mD3  u  D2  g  f4D4-  Let  (b,XQ,X1,p ,r2) 

((b' ,Xq,x|,p’ ,t*  ))  denote  the  parameters  associated  with  D  (D').  Then 

we  have  b  =  12m  +  15f4  +  6,  XQ  =  4m  +  5f4  +  2,  X1  =  f4  ^ 

b'  =  11m  +  15f4  +7,  Xq  =  3m  +  5f4  +  1,  X^  =  m  +  f4  +  1.  Thus  b  >  b' 
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when  m  ^  1.  We  have  p  <*  p'  iff  f4/(4m  +  6f^  +  2)  £  (m  +  f^  +  l)/(4m 

2  j  2 

which  holds  with  strict  inequality  when  m  >.  1.  Finally,  t  _>  t  iff 

4m  +  6 f  +2  4m  +  6f  +  2 

(4m  +  5f4  +  2)(4m  +  11  f^  +  2)  »  (3m  +  5f4  +  l)(9m  +  llf^  +  7) 

which  holds  with  strict  inequality.  This  completes  the  proof  for  Case  2 
hence  of  the  desired  result. 


< 


t  6f4  t  2) 


and 
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